Internal degrees of freedom and transport of benzene on graphite 
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In this paper, the chaotic internal degrees of freedom of a benzene molecule adsorbed on a graphite 
substrate, their interplay with thermal noise, and their effects on the diffusion and drift are investi- 
gated analytically by making use of the presence of two different time scales as well as by molecular- 
dynamics simulations. The effects of thermal noise are investigated, and it is found that noise does 
not significantly alter the dynamics of the internal degrees of freedom, yet affects the friction and 
diffusion of the center of mass. Qualitative and quantitative theoretical predictions for the friction 
and diffusion of the molecule on the substrate are made and are compared to molecular-dynamics 
simulations. Contributions to the friction and diffusion from the finite heat bath as well as the slow 
dynamics of the center of mass are formally identified. It is shown that the torsion in benzene, which 
dominates the nonlinear coupling, significantly affects the friction of the molecule on the surface. 
The results compare favorably with recent results from He/ncutron spin echo experiments on this 
system. Based on the analytical and numerical results, some suggestions are made for experimental 
conditions under which the effects of internal degrees of freedom might be observable. 
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I. INTRODUCTION 

The presence of mternal degrees of freedom (d.o.f.), es- 
pecially rotation, has been connected experimentally and 
theoretically to the diffusion and friction of nanoscale 
objects on surfaces (see, for instance, Refs. [IHS])- Re- 
cent He/neutron spin echo experiments [5HSj of benzene 
molecules adsorbed on graphite have shown that the fric- 
tion is surprisingly high and it has been conjectured that 
this could be related to the internal d.o.f. of the molecule. 

Several possible links can be identified between inter- 
nal d.o.f. and transport. The internal configuration of 
the molecules affects the interaction with the substrate, 
and therefore transport (see, for instance, Refs. p,llO)V 
However, the deterministic internal d.o.f. can also act as 
a source of noise [SJ [TT] , if their dynamics are sufficiently 
fast and chaotic. In this paper, I investigate the chaotic 
internal d.o.f. of a benzene molecule on a graphite sub- 
strate, their interplay with thermal noise, and their ef- 
fects on the diffusion and friction, by making use of the 
presence of two different time scales. The fast chaotic in- 
ternal degrees of freedom act as a finite heat bath which 
couples to the slow motion of the center of mass. I show 
that under realistic conditions, the time-scale separation 
theory of Ref. [12] is valid and links the dynamics of 
the internal degrees of freedom directly and quantita- 
tively to the diffusion and friction of the molecule on the 
substrate. The theoretical results compare well, quali- 
tatively and quantitatively, with simulations and exper- 
iments. One of the goals of this research is to construct 
systems where the friction or diffusion can be controlled 
by external means, by manipulating the internal degrees 
of freedom. 

In Sec. [nj the concept of time-scale separation in the 
context of nonlinear dynamics and surface adsorbates is 
introduced briefly. In Sec. |III[ the atomistic model for 
benzene on graphite used in this work is described. In 



Sec. |IV[ the time-scale separation theory is applied to 
the benzene on graphite model with thermal noise and 
expressions for the effective interactions are derived. It 
is also shown that thermal noise affects the dynamics of 
the internal d.o.f. weakly, but can still significantly affect 
the effective momentum diffusion and friction. A relation 
between the effective momentum diffusion and friction is 
also derived. Molecular-dynamics simulations of the sys- 
tem are described in Sec. Ivj in which specific d.o.f. can 
be frozen to investigate their effects, and thermal noise 
can be introduced slowly. The theoretical expressions are 
evaluated numerically and shown to compare favorably to 
simulation results and experimental results of benzene on 
graphite in Sec. |Vl] Finally, in Sec. |VII[ conclusions are 
drawn, and suggestions are made for experimental con- 
ditions under which the effects of chaotic internal d.o.f. 
might be observable. 



II. TIME-SCALE SEPARATION 

We first consider the nature of the coupling between 
the internal d.o.f. and the center of mass. In Refs. [Sllllj. 
the connections between transport and chaos in internal 
d.o.f. were outlined by taking advantage of the presence 
of multiple time scales. 

When a molecule like benzene is adsorbed weakly 
(physisorbed) on a surface like graphite, it is not strongly 
distorted by the surface. Accordingly, the internal forces 
are stronger than the forces exerted by the substrate and 
the time scales of the dynamics of the internal d.o.f. (for 
example vibrational modes) are much shorter than the 
time scales of the motion of the center of mass on the 
substrate. Specifically, the equations of motion can be 
split into those for the slow motion of the center of mass 
a = (R, P) and those for the fast internal coordinates 



(q,p) 



a = f(a,b) , 

b= g(a,b) 
e 



(1) 
(2) 



where f (a, b) and g(a, b) are both well-behaved functions 
of the same order of magnitude, and e is the time-scale 
separation parameter (0 < e ^ 1). 

To simplify the description of the slow subsystem, one 
usually wishes to eliminate the fast variables. Several 
strategies exist for doing this, depending on the dy- 
namical properties of the fast degrees of freedom. A 
(quasi)periodic fast subsystem can be averaged out [T5] 
while a highly disordered fast subsystem with many d.o.f. 
can be described as an infinite heat bath [T3]. 

The internal dynamics of complicated molecules, how- 
ever, are neither many nor (quasi)periodic 5j, but rather 
they are chaotic. A fast chaotic subsystem acts as a finite 
heat bath, leading to stochastic driving and damping of 
the slow system [T^ [T5HT7] (in this case the center of 
mass of the molecule on the substrate) , with clear signa- 
tures of the finite size of the fast system, and the finite 
amount of energy stored in it. 

When a dynamical system explores the entire phase 
space for almost all initial conditions, it is called ergodic. 
When, in addition, sets of initial conditions are smeared 
out after a long time, it is referred to as mixing. (See, 
for instance Refs. J18l[l9j.) If the fast subsystem has this 
property, its long-time behavior is independent of the 
initial conditions. If correlation in the fast subsystem 
also decays exponentially, then the dynamics of the slow 
coordinates a = (R, P) can be described by a Fokker- 
Planck equation for their probability density p(a, i) p^ . 



d 



dt 



5]^[4i)(a)p(a,0], (3) 



da^ 



where D'^^\&) is the 6-dimensional effective drift vector 
and D^'^\a) is the 6x6 effective momentum diffusion 
matrix. The subscripts /z, v indicate the six components 
of the slow coordinates, 1, 2, 3 referring to the positions 
and 4, 5, 6 to the momenta. The effective momentum 
diffusion and drift functions can be expressed in terms of 
averages over the invariant density of the fast subsystem 
for fixed slow coordinates, i.e. the probability density 
in the phase space of the fast subsystem which is left 
unchanged by the time evolution of the dynamics. One 
finds il2| 






dt'(<5/^(a,C(t',b;a))5/,(a,b))j 



0(e) 



(4) 
(5) 



DW(a) = (/4a,b))E 

(6) 



dttu 



= (/^(a,b))E + 0(e'), 



(7) 



where (.)e is used to indicate an ensemble average over 
the invariant density, ^(i, b; a) is the solution of the equa- 
tions of motion for the fast subsystem for fixed a and 
initial condition b, and 



5f(a,b)-f(a,b)-(f(a,b))E 



(8) 



The symbol O is used here to indicate the order of mag- 
nitude of the next terms, which are not necessarily neg- 
ligible. It should be noted that typical low-dimensional 
Hamiltonian systems are not mixing but have a com- 
plicated phase space, with regions of chaotic and quasi- 
periodic behavior |18H21| . However, the mixing prop- 
erty is not strictly required for the expressions above jllj . 
Furthermore, in Sec. |IV A[ it is shown that the addition 
of the substrate temperature ensures both exponential 
decay of correlation in the fast subsystem and mixing. 

In short, if the internal d.o.f. of a molecule are chaotic, 
then they act as a finite heat bath, providing noise to 
drive diffusion in the momentum [Eq. (HI)] and viscous 
friction [the second term in Eq. (l6|], affecting the diffu- 
sion and friction of the center of mass of the molecule on 
a surface |TT] . 



III. CHAOS IN MOLECULES: BENZENE ON 
GRAPHITE 

This work is concerned with the time-scale separation 
in benzene molecules on graphite, a prototype system for 
molecular adsorption. In order to investigate benzene, 
a model was developed in Ref. ^ which is particularly 
suitable for investigation of dynamical properties. 

Though molecules are in principle quantum- 
mechanical, often they can be described classically 
by atomistic models. Most typical atomistic force fields 
used for this (for example Tripos 5.2 [22) contain 
predominantly harmonic terms. Nevertheless, compli- 
cated geometries as well as torsion, both present in the 
internal degrees of freedom of benzene, lead to strong 
nonlinearities in the interaction and to chaos [Sj. 

The total potential energy of a benzene molecule con- 
sists of contributions from bending, stretching, and tor- 
sion of various types of bonds. However, due to the low 
mass of the hydrogen atoms, the energies associated with 
the vibrational modes involving the C-H bonds are high 
compared to the typical energy available at room tem- 
perature. Consequently, the C-H bonds will always be in 
the ground state, and do not participate in the chaotic 
dynamics O [23] . If an atomistic model is to be used to 
investigate the internal degrees of freedom of benzene, 
the C-H bonds must therefore be eliminated. Here, we 




FIG. 1. A schematic representation of a top view of a ben- 
zene molecule on a graphite substrate. Because the C-H bonds 
must be described quantum mechanically, they must be elimi- 
nated from the equations of motion before a classical atomistic 
description can be used. In the Tripos 5.2 force field three in- 
teraction parameters remain, the C-C stretching parameter 
kr, the C-CH-C bending parameter fc,^, and the C-CH-CH-C 
torsion parameter kp. 



average them out using a mean-field approximation, to 
obtain a greatly simplified Hamiltonian [5]. 

Let Ti and p; denote the position and momentum of 
the z-th CH complex, ordered in such a way that i and 
i + 1 are neighbors in the benzene ring. Let (j)i be the 
angle between r^-i — r.^ and r^+i — r^, and /3i the torsion 
angle associated with the bonds between the (« — l)st, 
ith, {i + l)st, and (i + 2)nd carbon atoms. The potential 
energy of the benzene molecule can be written as 



''benzene 1,^1 ; ■ • ■ ; ^6 J 



ifc,y^(||r(,+i)( 
i=l 


mod 6) 


6 

4=1 




6 

+ kp)_^[l + cos{2f3,)] , 



ro? 



(9) 



i=l 



where k^ and rg are the C-C stretching force constant and 
equilibrium distance, while k^ and kp are the effective 
bending force constant and the effective torsion constant. 
The model system is shown in Fig. [T] In the mean-field 
approximation of Ref. [S] for the Tripos 5.2 force field |22j, 
which will be used in this work, we have vq = 1.47 A, 
kr = 60.7 eV/A2, k^ = 6.85 eV/rad^ and kp = 0.247 eV. 
The typical energies (frequencies and thus time scales) 
can be estimated from the linearized system and range 
between 0.048 eV and 0.24 eV, which is comparable to 
the fundamental vibrational frequencies of benzene. 

Chaos can be quantified by the largest Lyapunov ex- 
ponent, the rate at which infinitesimal perturbations ex- 
pand. If this is positive, then the system is chaotic. The 
largest Lyapunov exponent of this model with energies 
equivalent to room temperature can be calculated from 



molecular-dynamics simulations and is A = 0.5/ps, re- 
gardless of the precise initial conditions [5] . The system 
is strongly chaotic and also (close to) ergodic and mixing. 



A. Graphite substrate with temperature 

Here, I consider benzene molecules adsorbed on a 
graphite substrate. Graphite is frequently chosen as 
a substrate in experiments. It is often modelled with 
two-dimensional sinusoidal potentials (see, for instance, 
Ref. [Tj). In this work a three-dimensional substrate po- 
tential is needed, and therefore the two-dimensional po- 
tential of a hexagonal graphite substrate in the xy plane 
is extended by adding a harmonic term in the z direction. 
The coefficient for the harmonic term is chosen such that 
in a potential minimum the second derivative of the po- 
tential energy is the same in all three directions. The po- 
tential to which the CH complex at position r — (x, y, z) 
is subject due to the substrate is given by Ij 



FcH(r) = - 



2V, 
9 



2 cos 



-t- cos 



27rx 
a\/3 



/2ny\ 



\ 3a J 



+ K 



27a2' 



(10) 



where V, = 25 meV is the potential difference associated 
with the corrugation, and a = 1.42 A is the in- layer inter- 
atomic distance of graphite. The origin and direction of 
the X and y axes are indicated in Fig. [TJ Note that the 
typical frequency of the dynamics of the center of mass 
in this potential (associated with the lowest energy mode 
at 0.38 meV) is more than two orders of magnitude lower 
than the lowest vibrational frequency of the molecule. 
The Hamiltonian of the system can now be written as 



CH 



1 ^ 

H = ^^^^^ ^ P^ + Vbcnzene(ri, . . . , Tg) 

6 

+ E^CH(r.), (11) 



where Pi is the momentum of the i-th CH complex, and 
men is the mass of a CH complex. 

In preliminary simulations without substrate thermal 
noise |5], it was shown that in this system the chaos in 
the internal d.o.f. does indeed lead to diffusive behavior 
of the center of mass and that this diffusion is normal, as 
expected from the time-scale separation. 

Because the substrate has a finite temperature, the 
molecule adsorbed on it is subjected to thermal fiuctua- 
tions. These are modelled here with Langevin dynamics 
applied to each CH complex: stochastic noise, driving 
the momentum with diffusion constant Dt, and viscous 
friction with momentum damping constant 77^, 



FLangevin = -iTmcHYi + rii{t) 



(12) 



where the amphtude of the noise must be chosen such 
that the expectation value of J^ dtr]^{0)-r]^{t) is equal to 
3Dt- The friction and diffusion constants can be related 
to the temperature T through 



Dt 

It 



= -mcakBT 



(13) 



The amplitude of the noise can be varied without chang- 
ing the temperature, by increasing 77- proportionally 
with Dt. Typically, the friction coefficient 7t of a 
graphite substrate is expected to be around 1/ps. How- 
ever, in recent He/neutron spin echo experiments, it was 
found that the friction coefficient of benzene on graphite 
was unexpectedly high, 2.2 ± 0.1/ps j5j. 



IV. EFFECTIVE INTERACTION 

The effective dynamics of the benzene molecule on the 
substrate are determined by the coupling and noise from 
the substrate and molecule. In this section, the effective 
interaction is investigated analytically in more detail, and 
the various contributions are formally identified. 

In order to apply the time-scale separation, and use 
the expressions of Eqs. Q and ^ to determine the ef- 
fects on the diffusion and friction of the molecule, the 
equations of motion must first be rewritten in the form 
of Eqs. ^ and ([2|, by separating out the slow coordi- 
nates R, P of the center of mass from the 15 fast internal 
d.o.f., which are denoted by qi, . . . , qs, with associated 
momenta pi, . . . , ps. We have P = 6mcHR- 

Let the lowest frequency of the linearized benzene 
molecule be denoted by Wbenzenej and the typical fre- 
quency of the substrate be denoted by Wgraphitc- The 
time-scale separation parameter is related to these two 
frequencies (given in the previous section) via 



e = """"P^'*^ « 0.0078 

'^benzene 



(14) 



Let the set of internal coordinates be denoted by 
r, R 



V~e ' 



(15) 



for z = 1, . . . , 5. The corresponding momentum is 

Pj = x/emcH(i"j - R) • (16) 



For convenience of notation, we shall use the notation 
qe = -SjLiQi and pe = -^LiPj- ^^ the potential 
energy of benzene, Eq. (lo]), contains only differences of 
the positions of the CH complexes, and not R explicitly. 



Eq. ( 11 ) can be written as 



12mcH 2toch e 



1 



1=0 



+ Vbcnzonc(qi/%/e, • • ■ , Qe/Ve) 

6 

+ ^ycH(R + q./^) , (17) 

i=l 

The equations of motion of the full system read, 

R=^r^. (18) 



6mcH 
d 



P = - aS E ^Ch(R + q./^) -1TP + J2 ^^(*) ' 
i=l i=0 



2mcHe 
d 

Pi = - T^- ^benzene (qi/\/e, • • • , (Ig/V^) 



(19) 
(20) 



9q.i 



_d_ 
9q, 



FcH(R + q»M) 



iTPt + Ve \ Vt{t) - iYl ^j(^) 



(21) 



From a comparison between the energy constants of the 
benzene molecule [Eq. (IqI)], and the potential of the sub- 



strate [Eq. ( 10 1] , which is at least two orders of magnitude 
smaller, it can be concluded that the right-hand sides of 
Eqs. (20) and (21) are at least two orders of magnitude 



larger than the right-hand side of Eq. ( 18 ) and the first 
term on the right-hand side of Eq. (19). Consequently, 



as long as the Langevin noise and friction due to the sub- 
strate temperature have the same slow time scales as the 



rest of the substrate, Eqs. (18) and (19) are of the form 
of Eq. ^, while Eqs. ([20| and ^ are of the form of 
Eq. (|2|. 

Equations ( 18 21 ) can now be combined with Eqs. ([4| 
IgI) to calculate the effective interaction. Let (5F(R;q, 
be the deviation of the total coupling F(R; q^) from the 
average, 

SF{R; q,) = F(R; q,) - (F(R; q,0)E , (22) 



F(R; q. 



-^^ycH(R + q.M). (23) 



Equation ( 18 ) does not contain the fast variables, and 



therefore only the momentum of the center of mass is 
directly affected by the noise from the fast degrees of 
freedom. This means that only components of Dj^J with 
fi,iy > 3 can be nonzero. The diffusion due to the noise 
from the two heat baths can then be written as 



^(2) ^ j £'^z.(R) + Sfj,i,DT ii n,iy > 3 
^"^ I otherwise 



^p.(R) 



dt{SFiR;a^^,t))SF{R;q,))E 



(24) 
(25) 



where -D(R) is the effective momentum diffusion matrix 
due to the finite heat bath of the internal d.o.f. and 5^^ is 
the Kronecker delta. The solution of equations of motion 
for the fast coordinates with initial conditions q^ and 
fixed slow coordinates is denoted by ({q_i,t). 



^H^ if M < 3 

-Ea^T^^^^^ + '^'/^V - ^A' • 7(R) - 7tPm ifM>3 



(26) 



where 7(R) is the 3x3 effective friction matrix, which 
is strictly positive definite. It can be determined from 
the couphng and the decay of correlation in the internal 
d.o.f. by 



it.A'R) 



1 



6mcH 



dt 



JF^(R;C(q.,0) 
9R. 



(27) 



The first term in the second line on the right-hand 
side of Eq. (26), the average coupling due to the sub- 



strate potential, can be thought of as a configurational 
term. It contains the invariant density of the internal 
d.o.f. through the ensemble average, which can be af- 
fected in detail by the dynamics of the internal degrees 
of freedom. The second term is the viscous friction term 
due to the decay of correlation in the finite heat bath, 
while the third term is the viscous friction due to the 
coupling to the infinite heat bath of the substrate. Both 
the friction and diffusion depend on the location on the 
substrate and the direction of the momentum. 

Eqs. (25 1 - (27) are a central result, giving theoreti- 



cal expressions for the reduced dynamics of the center of 
mass. In Sec. |VI[ these expressions will be evaluated and 
compared to simulations. 



A. The effects of thermal noise on the internal 
d.o.f. 

The thermal noise from the substrate temperature not 
only affects the center of mass motion, but also the fast 
dynamics in the internal degrees of freedom. As it is 
uncorrelated Gaussian white noise, it ensures exponential 
decay of correlation in the fast subsystem. 

Because the coupling to the substrate is weak com- 
pared to the coupling between internal d.o.f., the fast 
subsystem is essentially subject to very weak Gaussian 
white noise and weak damping [see Eq. (21)]. Conse- 



quently, the dynamics remain similar to the dynamics 
without noise. 

Nevertheless, the noise still has important qualitative 
and qualitative effects on the invariant density, and thus 
indirectly on the ensemble averages and the effective mo- 
mentum diffusion and friction. In particular, the noise 
smears out the dynamics and structures in the phase 
space [20]. This is important, because the results of 
Ref. [12] require that the fast subsystem be mixing, 
while low-dimensional Hamiltonian systems are, in gen- 
eral, not. However, with a suitable realization of the 



Gaussian noise, any point in phase space can be reached 
in finite time from any initial condition. In other words, 
the entire phase space is explored in finite time. Together 
with the smooth dynamics of the benzene molecules, the 
above guarantees mixing. Consequently, a unique invari- 
ant density exists, and the conditions of the time-scale 
separation theory are met. 

In the presence of thermal noise, the fast subsystem 
is simply coupled weakly to an infinite heat bath, with 
a temperature equal to the substrate temperature T. In 
the absence of other external driving, the ensemble av- 
erage becomes equal to the thermal average, with the 
invariant density equal to the Gibbs distribution, 



(A(R,P,qi,...,p5))E 



— / rfqi . . . (ip5 exp 



H 

kWT 



A(R,P,qi,...,p5) . 



(28) 



Weak noise can change the invariant density significantly, 
and therefore the reduced slow dynamics through the av- 
erage coupling, as well as the effective friction and diffu- 
sion. 

In addition to this, it is possible to obtain a direct 
relation between the diffusion and drift. An addition 
coupling term can be applied to the center of mass, which 
restricts it to a small region around a particular position 
R on the substrate. The internal dynamics will not be 
affected, and only the average interaction [the first term 
in the second line of Eq. (26)] term will change. The 



effective friction and diffusion will therefore remain the 
same, and can be approximated as constants if the region 
is sufficiently small. Hence, the center of mass is simply 
coupled to an infinite heat bath which produces diffusion 



and friction tensors D^,^ and G^j^ with 



-•HD 



^^,(R) + i^T<5p. , 

= 7mi.(R) +lT^t,v ■ 



(29) 
(30) 



The center of mass is coupled to an infinite heat bath 
consisting of the substrate and internal degrees of free- 
dom of the molecule. This heat bath necessarily has the 
same temperature as the substrate. Consequently, sim- 



ilar to Eq. ( 13 ) but in three dimensions, and allowing 
for the possibility of negative eigenvalues of the diffusion 
tensor, 



7a.. (R) 



i^Si^" 



;R)|d«^(R)d,,(R), (31) 



where I?k(R.) is the Kth eigenvalue of I?^y(R), and 
d^t (R) is the bth component of the corresponding eigen- 
vector. Equation (31) is a particularly useful result, as 
Eqs. (25) and (27) contain different information regard- 



The equations of motion based on Eq.(ll), with 
Eqs. ([9]) and (10) substituted, were solved using the 



ing the symmetries of the diffusion and friction tensors. 
These symmetries can be used to improve the numerical 
evaluation of the theoretical expressions, and are used for 
this purpose in Sec. |VI[ 



B. Transport on the substrate 

The diffusion and friction tensors cannot, as yet, be 
measured directly in an experiment. What we are inter- 
ested in, therefore, is the combined average effects of all 
contributions to friction and diffusion, i.e. decay of cor- 
relation, on a typical trajectory on the substrate. Besides 
the loss of correlation due to the thermal noise and the fi- 
nite heat bath of the internal d.o.f. of the molecule, corre- 
lation can also decay due to chaotic dynamics in the slow 
coordinates. As both diffusion coefficients and friction 
coefficients are additive, and the substrate is symmetric 
under rotations of 60°, this can be written as 

7total,;ii/ — <5^iy(7slow + Tbcnzcne + 7t) , (32) 

Dtotal,fj.i' ~ <5/jiy (-Dslow + ^'benzene + Dt) ■ (33) 

where /z, i' > 3, and the subscript slow and benzene indi- 
cate the contributions from the slow coordinates and in- 
ternal d.o.f. of the molecule respectively. Note that if the 
substrate lattice were rectangular, instead of hexagonal, 
the average friction and diffusion would not be isotropic. 
As the invariant density is not affected by the noise 
level, and the fast dynamics are not affected strongly for 
up to reahstic noise levels (1/ps), 7bcnzcno and L'bonzcno 
are constant for fixed temperature, irrespective of jt- 
The slow dynamics are affected directly by the thermal 
noise, but this effect is small up to realistic noise levels. 
They are also affected indirectly by the presence of the 
thermal noise, but not its intensity, through the invariant 
density, as can be seen from the first term in the second 
line of the right-hand side of Eq. (26), which contains 



the ensemble average. The friction and diffusion of the 
reduced slow system, 7siow and I?siow are not straight- 
forward to treat analytically, even though this system 
is essentially a single particle moving in a complicated 
substrate potential (see, for instance, Ref. [M])- 



velocity- Verlet algorithm or in the case of Langevin 
damping and noise with a fourth order Runge-Kutta 
algorithm. As this work is concerned with the dy- 
namical properties of the system, it is necessary in the 
Hamiltonian case to use an algorithm which preserves 
the phase space volume exactly. The largest Lyapunov 
exponent was obtained by numerically calculating the 
growth rate of an infinitesimal perturbation. By virtue 
of the Oseledec theorem, this growth rate is equal to the 
largest Lyapunov exponent (see, for instance. Chap. 3 of 
Ref. ES]). 

Random initial conditions for the center of mass and 
the internal d.o.f. were chosen with energies correspond- 
ing to room temperature. This was done by starting from 
a random point on the constant energy shell of the lin- 
earized system, with energy corresponding to room tem- 
perature and zero angular momentum, and subjecting 
the system to thermal noise for a period of time suffi- 
cient for reaching thermal equilibrium. 

The total friction coefficient was extracted from the 
trajectories by fitting a linear dependence to the aver- 
age change in velocity as a function of the momentum, 
(P(t+At)— P(i))traj- Note that in order to accurately de- 
termine the friction due to the internal d.o.f.. At must be 
sufficiently long for the correlation in the internal d.o.f. 
to have decayed. Similarly, the diffusion of the center 
of mass of the molecule on the substrate can be deter- 
mined from the mean square displacement ([R(t -|- At) — 
R(t)]^)traj as a function of the length At of the time in- 
terval, as was done in, for instance, Refs. [11] [55]. If the 
diffusion is normal, ([R(i + At) — R(i)]^)traj is propor- 
tional to At for long At, and the prefactor is proportional 
to the diffusion coefficient. If the mean square displace- 
ment does not grow linearly with time, but with some 
other exponent a, 



{[Kit + At) - Kit)] 



} traj 



(Xt° 



(34) 



the diffusion is called anomalous. Anomalous diffusion 
can be subdiffusion (a < f) or superdiffusion [a > 1). 
Note that anomalous diffusion cannot be described by a 
linear Fokker-Planck equation such as Eq. (Isl). 



Fixed internal d.o.f. 



V. MOLECULAR DYNAMICS SIMULATIONS 

Molecular dynamics simulations were performed of the 
model system, and total friction and diffusion constants 
were calculated numerically. In the simulated system, 
specific internal d.o.f. can be frozen, and the thermal 
noise due to the substrate temperature can be introduced 
slowly, thus providing insight into the effects of the in- 
ternal d.o.f. and substrate temperature on transport. 



Besides the full system, several derived systems were 
considered with different sets of frozen internal d.o.f. with 
progressively fewer remaining internal d.o.f. of the ben- 
zene molecule, namely (1) frozen torsion modes, (2) all 
internal vibrations frozen while leaving rotation, and (3) 
fully frozen internal d.o.f. with (a) the molecule at in- 
commensurate orientation, (b) the molecule at commen- 
surate orientation, and (c) thermal average over all orien- 
tations. Commensurate orientation means that the ben- 
zene molecule has the same orientation as the hexagons of 



TABLE I. Several dynamical properties for the full system 
and various sets of frozen internal d.o.f., namely the system 
without torsion, without any vibration, and without any in- 
ternal d.o.f., for a fully frozen molecule at commensurate and 
incommensurate orientations, and for a thermal average over 
the orientation. Shown are the largest Lyapunov exponents A 
of the noiseless system with and without a substrate in units 
of 1/ps at energies comparable to room temperature, the ex- 
ponent of the diffusion q, and the induced friction 7 — 7t 
when the Langevin dynamics are switched on. 
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frozen, average 


(3c) 
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< 0.6 


1 < a 


< 2 


0.74 



C atoms in the substrate. The results which are described 
in the rest of this section are summarized in Table ID 

The lowest three fundamental vibrational frequencies 
of the benzene molecule correspond to torsion modes, 
in which the C atoms oscillate in the z direction. The 
other modes have higher frequencies and involve exclu- 
sively bending and stretching. By removing the dynam- 
ics in the z direction, the torsion modes can be elimi- 
nated while leaving the other vibrational modes intact 
[case (1)]. Removing the torsion modes greatly affects 
the chaos in the internal degrees of freedom, because 
their coupling is strongly nonlinear and dominates the 
chaotic dynamics. It also decreases the time-scale of the 
internal d.o.f., thereby increasing the separation between 
the time scales and decreasing e. In molecular dynam- 
ics simulations without temperature, with a fixed center 
of mass position, and without the torsion modes, it was 
found that the largest Lyapunov exponent is typically 
below 0.014/ps depending on initial conditions. This is 
much smaller than the value for the full system with fixed 
center of mass, 0.5/ps. Without torsion, the largest Lya- 
punov exponent is, for some initial conditions, even zero 
within numerical error. The resulting trajectories are 
not appreciably chaotic. Note that due to the small time 
steps used in the simulations, the error in the Lyapunov 
exponents for a particular set of initial conditions is less 
than 0.001/ps. There is, however, variation due to the 
variation in the initial conditions. 

The bending and stretching modes have higher fre- 
quencies and weaker nonlinear coupling. These modes 
can be frozen as well, by making the molecule completely 
rigid, leaving only one internal d.o.f., the rotation around 
an axis parallel to the z axis [case (2)]. With only one 
internal degree of freedom, for fixed center of mass coor- 
dinates the molecule cannot be chaotic. Because the fast, 
strong coupling between the atoms has been removed, the 
remaining dynamical system no longer has a fast time 
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FIG. 2. Sinrulated trajectories of the center of mass of a ben- 
zene molecule on graphite without temperature, for the full 
system, as well as various systems in which some of the d.o.f. 
have been frozen. In all cases the simulations lasted for 72 ns 
and initially the center of mass was at position R = (0, 0, 0). 
Two distinct trajectories are shown for the full system, two 
trajectories for a system without motion along the z direction, 
and two trajectories for a system with a rigid molecule which 
is allowed to rotate. For fully frozen benzene molecules one 
trajectory is shown with incommensurate orientation and one 
with commensurate. Long ballistic jumps can be seen even 
for the full system, which has normal diffusion. 



scale. The time scale of the rotation is comparable with 
that of the center of mass motion. Note that in very large 
rigid nanoclusters, the rotation is much slower than the 
center-of-mass motion [31 [^^ . 

For fully frozen internal d.o.f., the molecule is rigid 
[case (3)]. Depending on the initial orientation, the in- 
teraction with the substrate is different. Therefore, both 
commensurate and incommensurate initial orientations 
are considered, as well as systems where the interaction 
is equal to the thermal average for a rigid molecule that 
can still rotate. This is essentially equal to the second 



line of Eq. ( 26 1 with the expression for the ensemble av- 



erage in noisy systems, Eq. ( 28 ) for a rigid molecule with 



only rotational degrees of freedom, and a good first ap- 
proximation of the thermal average for a molecule with 
all internal degrees of freedom. 



B. Trajectories 

In Fig.[2j examples of simulated trajectories of the cen- 
ter of mass on the substrate of equal duration are plotted 
for different dynamics of the internal d.o.f. and without 
thermal noise. The full system is chaotic, with the largest 
Lyapunov exponent typically around 1.2/ps regardless of 



the precise initial conditions. The motion of the center 
of mass on the substrate is diffusive, but is dominated 
by intermittent long-range ballistic motion. The diffu- 
sion of the center of mass is normal with a total diffusion 
constant of around 27 nni^/ns. 

When the torsion modes are eliminated [case (1)], or 
when the internal d.o.f. are frozen completely except 
for rotation [case (2)], the dynamics remain chaotic. 
The largest Lyapunov exponent in both cases is around 
1.2/ps, similar to the full system. The motion of the cen- 
ter of mass is diffusive. For many initial conditions the 
diffusion is normal on the time scales of the simulations 
(120 ns), but for some initial conditions, superdiffusion 
cannot be excluded. As the molecule itself in these two 
systems is either not chaotic, or very weakly chaotic, it 
is the addition of the substrate that leads to chaos. In 
low-dimensional Hamiltonian systems with such a com- 
plicated phase-space structure, anomalous diffusion is to 
be expected for some initial conditions (see, for instance, 
Ref . [23 and some of the chapters in Ref . [55] ) ■ 

When all internal d.o.f. are removed [case (3)], the 
system becomes more weakly chaotic, with the largest 
Lyapunov exponent below 0.6/ps. For some initial con- 
ditions, there is still diffusive behavior. In nearly all of 
these cases, the diffusion is clearly superdiffusion, with 



scahng exponents a up to 1.9 [see Eq. (34|]. 

In the system under consideration here, the diffusion 
is not governed by single hops, but by ballistic motion in 
the directions that correspond to low diffusion barriers, 
as can be seen from Fig. [2j Ballistic motion was also 
observed in experiments [B]. The anomalous diffusion 
in the systems without torsion is also associated with 
these ballistic jumps. In the full system, the noise from 
the chaotic dynamics of the torsion modes causes the 
ballistic trajectories to decay more quickly and reduces 
the total diffusion of the center of mass on the substrate 
to normal diffusion. However, it is still dominated by 
long jumps, and, on short time scales, the diffusion is 
still anomalous, as has for instance also been found for 
single particles with Langevin dynamics in periodic and 
random potentials [26] [29] . 

The results for the Lyapunov exponents and diffusion 
exponents are summarized in Table llj 



C. Decay of correlation in the internal d.o.f. 
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FIG. 3. The auto-correlation function for various compo- 
nents of the coupling between the center of mass and internal 
d.o.f. as a function of time at R = (0, 0,0),P = (0,0,0) for 
a single simulated trajectory. Due to the symmetries of the 
lattice in point (0,0,0), the auto-correlation functions for y-y 
coincides with that for x-x. From the theory, all elements 
must integrate to zero. Modes with lower frequencies decay 
faster. The inset shows a semi-log plot of the second deriva- 
tive of the x-x auto-correlation function, which emphasizes 
the high frequencies, and falls off exponentially. The dotted 
line in the inset is an exponential fit to the maxima with a 
decay time of 0.55 ps. 



In Fig. [3] the auto-correlation function of the coupling 
is plotted for a benzene molecule kept at fixed position 
R = (0, 0, 0) and momentum P — (0, 0, 0). It decays ex- 
ponentially with a decay time of 0.55 ps. The vibrations 
which persist for a long time correspond to the high- 
est vibrational frequencies of the linearized system. As 
the nonlinear coupling of these modes to other modes is 
weaker than that of the modes with lower frequency, cor- 
relations in them decay more slowly. The time-integral 
of the autocorrelation function gives an estimate of the 
velocity-diffusion coefficient, via Eq. Q, which is equal 
to 0.0009 A^/ps^. This value is close to 0, which is the 
value expected from Eq. (31 1 in combination with the 
^ (0,0,0). 



symmetries of Eq. (27) in R = (0,0,0). As the internal 
dynamics of the molecule are not strongly affected by 
the substrate, the decay in the internal d.o.f. is similarly 
exponential for all positions on the substrate. 



The effective drift and diffusion are affected by the cor- 
relation times of the internal d.o.f., as the effective mo- 
mentum diffusion and friction due to the internal d.o.f.. 



as expressed in Eqs. (25) and (27), contain integrals 
over correlation functions of the coupling. The auto- 
correlation function in the expression for the momentum 
diffusion, the integrand on the right-hand side of Eq. ( 25 ) , 



is straight-forward to evaluate numerically. If the inte- 
gral over the auto-correlation function is not finite, the 
diffusion, if any, is anomalous and cannot be described 
by a linear Fokker-Planck equation such as Eq. (|3|. 



D. Thermal noise 

We introduce thermal noise and vary the intensity 
between zero to realistic values, while setting the tem- 
perature equal to room temperature, as described in 
Sec. |III[ In Fig. |4j the diffusion coefficient of a ben- 
zene molecule with various types of internal d.o.f. on 
the graphite substrate at room temperature is shown 
as a function of the imposed thermal friction 77^. As 
soon as the thermal noise is switched on, correlations 
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FIG. 4. The diffusion coefficient of the center of mass on the 
substrate is plotted as a function of the thermal noise inten- 
sity (characterized by the imposed friction 7t) in molecular- 
dynamics with Langevin dynamics at a temperature of 293 K. 
Without Langevin dynamics, the systems which are totally 
frozen do not exhibit normal diffusion. The error in the ef- 
fective diffusion in the reduced system at low noise intensities 
is large, around factor of 3. The diffusion coefficient of the 
full system without thermal noise is equal to 27 nm'^/ns. The 
diffusion coefScients were obtained from trajectories with a 
total duration of 2 ps. 



in the fast subsystem decay exponentially, and diffu- 
sion becomes normal, regardless of the dynamics of the 
internal degrees of freedom. Long ballistic jumps are 
destroyed faster by higher noise intensities, leading to 
slower diffusion. At the realistic intensities correspond- 
ing to 7t = l/ps, the diffusion of the center of mass for 
all systems which have the rotational degree of freedom 
is around 6.9 nm^/ns, which, given the non-empirical na- 
ture of the substrate potential, is close to the experimen- 
tal result of 5.39 ± 0.13 nm^/ns [6]. 

With noise, the full system, the system without tor- 
sion, and the system without any vibration have very 
similar diffusion constants. In these systems, the diffu- 
sion is slower than in systems without any internal d.o.f., 
due to the faster destruction of ballistic jumps. In sys- 
tems in which the diffusion is not dominated by ballistic 
motion, but by single hops between minima in the sub- 
strate potential, the additional noise due to the chaos 



in the internal d.o.f. might have more non-trivial effects, 
potentially enhancing the diffusion of the center of mass 
compared to thermal noise alone. 

Of the three systems without internal d.o.f., the system 
with the molecule at incommensurate orientation has the 
highest diffusion coefficient, because it has the weakest 
interaction with the substrate. At jt ~ 0.25/ps, the dif- 
fusion in systems with internal d.o.f. and systems without 
internal d.o.L and average interaction has become nearly 
the same. At this point, the noise due to the chaotic 
dynamics becomes weak compared to the thermal noise, 
and no longer affects the diffusion constant significantly. 
As the expected noise intensity in the physical system 
corresponds to around l/ps, this means that the diffu- 
sion of the benzene molecule is not significantly affected 
by the noise due to the finite heat bath of the inter- 
nal degrees of freedom. However, in systems where the 
time-scale separation parameter e is larger, the effective 
momentum diffusion will be larger [see Eq. (HI)], and may 
become significant. Furthermore, in the experiments of 
ReL P, long jumps dominate the diffusion even at realis- 
tic noise intensities, and consequently noise from internal 
degrees of freedom can modify diffusion significantly. 

In Fig. [5| the induced friction coefficient is plotted as 
a function of the imposed friction. The time interval 
used to estimate the friction was 0.484 ps, as most of the 
correlations in the internal d.o.f. have decayed after this 
time (see Fig. [s] and Sec. |VC[ ). At jt = 2.0/ps, the total 
damping during the time interval of 0.484 ps is so strong 
that it becomes difficult to obtain an accurate estimate 
of the friction. However, as can be seen from Fig. [3j for 
shorter time intervals correlation has not decayed suffi- 
ciently to obtain accurate estimates. 



VI. COMPARISON BETWEEN THEORY, 
SIMULATIONS, AND EXPERIMENT 

The qualitative theoretical results for the relation be- 
tween internal chaos and diffusion compare well to the 
simulation results, as is summarized in Table |l] In the 
full system, where the internal d.o.f. are chaotic (strictly 
positive Lyapunov exponent), and correlation decays ex- 
ponentially, the noise from the internal d.o.f., combined 
with the time-scale separation, lead to normal diffusive 
behavior. When the internal d.o.f. arc not chaotic for 
some or all initial conditions, diffusion can be anoma- 
lous. In the rest of this section, also quantitative results 
from the theory for the reduced system are compared to 
results from simulations of the full system. 



Evaluating the theoretical expressions 
numerically 



The theoretical expressions given in Eqs. (25) and (27) 



contain averages and integrals over correlation functions 
of the fast subsystem for fixed slow coordinates. These 
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FIG. 5. The component of the friction due to the interaction 
between the internal d.o.f. and the substrate is plotted as a 
function of the imposed friction in molecular-dynamics simu- 
lations with Langevin dynamics at a temperature of 293 K. 
The friction coefficients were obtained by averaging over tra- 
jectories with a total duration of about 2 fis. For comparison, 
the friction constant of the full system without thermal noise, 
7o ~ 0.75/ps, is also plotted. 



expressions are difficult to evaluate in any more detail 
analytically, but here are evaluated numerically by per- 
forming simulations of the internal d.o.f. for fixed slow 
coordinates and calculating the necessary averages and 
correlation functions. In these numerical simulations, 
the parameter jx was set to a non-zero value, 0.1/ps, 
to ensure that the system is mixing and has the correct 
invariant density. 

The average coupling term [the first term of the second 



line of Eq. ( 26 )] is simple to evaluate in numerical simula- 
tions. Similar y, the auto-correlation function needed for 
the momentum diffusion, shown in Eq. ( 25 ) , is straight- 



forward to obtain from simulated trajectories and to in- 
tegrate. The results for the effective momentum diffusion 
as a function of the position are shown in Fig. [6] 

The expression for the friction coefficient [Eq. (27l] 
is more complicated to evaluate quantitatively, due to 
the derivative. Because the invariant measure depends 
on the position of the center of mass, the average and 
the derivative are not interchangeable. Hence, in order 
to calculate the derivative, the average coupling must 
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FIG. 6. The effective velocity diffusion coefficient D^^(R) 
as a function of position, obtained by evaluating the theoret- 
ical expression, Eq. (251, numerically. For reference, carbon 



atoms and bonds between them are indicated with solid circles 



and dotted lines. From Eq. (251 it is clear that the effective 



momentum diffusion does not depend on the velocity or on 
the z coordinate. 
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FIG. 7. The effective friction coefficient 7;j^(R) as a function 
of po siti on, obtained by evaluating the theoretical expression, 
Eq. (271, numerically. Autocorrelations were calculated for 
up to 0.7 ps and integrated. For reference, carbon atoms 
and bonds between them are indicated with solid circles and 
dotted lines. As with the diffusion, the effective friction coeffi- 
cient does not depend on the velocity or the z coordinate [see 
Eq. (27l]. As it is more computationally involved to calculate 



the friction tensor than the diffusion tensor, the results are 
less accurate. Within the error of the numerical evaluation, 
the effective friction indeed has the same shape as the diffu- 
sion, shown in Fig. [6] with the correct ratio, and thus obeys 
Eq. pll). 
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be evaluated sufRciently accurately at nearby positions. 
Numerical evaluation of this expression is therefore more 
time consuming than for the diffusion, and less accurate. 
In order to obtain sufficient statistics, the nearby posi- 
tions had to be chosen to be 0.26 A apart. Results for 
the effective friction coefficient are shown in Fig. [7| 

As can be seen from Figs. [6] and [7) within the error of 
the numerical evaluation, the effective momentum diffu- 
sion and friction indeed have the same shape, with the 



dominated by long ballistic jumps, as was mentioned in 



correct ratio, and thus obey Eq. (31). 



B. Comparison of theory with simulations of the 
full system 

It is possible to compare the theoretical results to 
numerical trajectories. However, comparing the results 



from Eqs. (25) and (27 1 that were obtained numerically 



in the previous subsection to diffusion and drift tensors 
obtained directly from numerical simulations of the full 
system requires an accurate determination of the mean 
square displacement in velocity space from trajectories of 
the entire system. However, due to the high dimension- 
ality of the slow subsystem, extremely long trajectories 
are needed to obtain useful statistics. With currently 
available computing power this is not yet possible. It 
is interesting to note that the theory can thus be used 
to calculate properties of the system that cannot yet be 
obtained from molecular-dynamics simulations. 

Nevertheless, the theory can be tested by comparing 
total diffusion and friction coefficients of the center of 
mass on the substrate found from numerical simulations 
of the reduced dynamics to those obtained from simula- 
tions of the full system. 

The numerical values for the theoretical expressions 
at specific points, obtained in the previous subsection, 
can be used to obtain an approximation for the effective 
drift and diffusion, by fitting it to a Fourier expansion 
with the appropriate symmetries, cut off after a certain 
wave vector (in this case Sir /a). These symmetries are 
the translation, rotation, and reflection symmetries of the 



lattice, the reflection symmetries of Ec^. (27), the symme- 



try of the diffusion tensor in Eq. ( 25 1 , and finally Eq. ( 31 1 



Though, in principle, there are many terms, due to the 
large number of symmetries only four free parameters 
remain. 

Results for the total friction and diffusion of the re- 
duced system are included in Figs. |4] and [5j At low noise 
levels, the induced friction in the full system is about 
0.83/ps, while the reduced system has an induced fric- 
tion of (0.86 ± 0.02)/ps. The error is due to the error 
in the (fit to the) numerical evaluations of the effective 
momentum diffusion and friction and was estimated with 
a variational method. 

The diffusion of the center of mass in the reduced sys- 
tem follows that of the full system for high noise levels. 
For low noise levels, however, there is a discrepancy. This 
is due to the fact that in this regime, the diffusion is 



Subsec. V B The long jumps are very sensitive to noise 
and damping, and disappear at high noise levels. The 
noise levels needed to destroy the long jumps are com- 
parable to the typical intensities of noise and damping 
from the internal degrees of freedom. Nevertheless, the 
long jumps survive in both the full and reduced system, 
because for some positions on the substrate there are no 
noise and damping from the internal degrees of freedom. 
The duration of the jumps, however, depends sensitively 
and in a highly nontrivial way on the reduced dynamics 
in these regions, where the noise level due to the internal 
d.o.f. is extremely low. Due to the statistical errors in 
the numerical evaluation of the theoretical expressions, 
the error in the effective momentum diffusion at these 
noise intensities is large, around a factor of 3. It should 
be noted that while the thermal damping jt is usually 
taken to be a constant, it too could, in principle, depend 
on position, which could affect the survival of long jumps. 



C. Various internal degrees of freedom 

As the intensity of the thermal noise affects the dy- 
namics of the internal d.o.f. only weakly, and not the 
invariant density, the effective friction and momentum 
diffusion are not affected by it. From Fig. [5] it can be 
seen that the induced friction, which contains the effec- 
tive friction due to the finite heat bath of the internal 
d.o.f., does not depend on the noise intensity. Together 
with the values for the induced friction in the various sys- 
tems, one can use this to separate out the contributions 
to from the internal degrees of freedom. 

The substantial difference between the induced fric- 
tions in the fully frozen systems with commensurate and 
incommensurate orientations (0.99/ps and 0.29/ps re- 
spectively) demonstrates that the interaction between 
the internal d.o.f. of the molecule and substrate strongly 
affects the friction. The full system and other systems 
with reduced or partially frozen dynamics lie somewhere 
in between these two extremes. The next best approxi- 
mation, the system without any internal d.o.f., but with 
the coupling averaged with thermal weight over orienta- 
tion has an induced friction coefficient of 0.74/ps. Be- 
cause it is not related to the internal d.o.f., this amount 
can be attributed to the slow dynamics. 

The induced friction in the full system is 0.83/ps for a 
wide range of noise intensities, whereas, without thermal 
noise the full system has an induced friction of 70 = 
0.77/ps. This difference is due to the fact that even weak 
thermal noise changes the invariant density significantly 
compare to the system without noise (see Sec. IV A) and 
this affects the friction induced by the finite heat bath 
[see Eq. ^]. 

The induced friction in the systems which still possess 
rotation, but have frozen torsion modes, is 0.77/ps in 
both the system with and without bending and stretch- 
ing. The difference between this value and the value for 
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the full system can be attributed to the chaotic dynam- 
ics of the torsion modes. The remaining difference to the 
system with a rigid molecule with thermal average over 
orientation is due to rotational dynamics. 

In summary, three contributions to the total friction 
can be identified: the damping due to the thermal heat 
bath of the substrate (77- w 1/ps), a constant contribu- 
tion from the slow dynamics (7siow = 0.74/ps), and the 
contribution from the chaotic internal d.o.f. (7benzcnc = 
0.09/ps), which is dominated by the torsion modes. This 
contribution from the internal d.o.f. is a combination of 
the average coupling and finite heat bath due to the in- 
ternal vibrations. 

The total induced friction of 0.83/ps found in simula- 
tions matches well with the theoretical predictions from 
the time-scale separation theory, 0.86±0.02/ps. With the 
addition of the 1/ps substrate friction due to dissipation 
into phonon modes, the numerical and theoretical calcu- 
lations compare favorably with the result of 2.2 ± 0.1/ps 
from Hc/ncutron spin echo experiments [B]. 



VII. CONCLUSIONS 

In this work, the relationship between the internal 
d.o.f. of a benzene molecule and its motion on a graphite 
substrate were investigated using time-scale separation 
theory and dynamical properties, including the effects 
of thermal fluctuations. It was shown that the fast dy- 
namics of the internal d.o.f. act as a flnite heat bath, 
while the motion of the center of mass is slow. It was 
found that thermal noise affects the internal dynamics of 
the molecule weakly, yet ensures the requirements of the 
time-scale separation theory of Ref. [12], i.e. exponential 
decay of correlation and mixing. The invariant density 
of the internal d.o.f. in the presence of noise is equal to 
the thermal distribution. The total friction and diffusion 
of the center of mass consist of contributions from the 
thermal heat bath of the substrate, the finite heat bath 
of the molecule, and the decay of correlation in the slow 
dynamics. 



The theoretical expressions for the effective dynamics 
were evaluated numerically, and numerical simulations of 
the reduced system were performed. Their results com- 
pare well with those from simulations of the full system. 
In the molecular dynamics simulations, several sets of in- 
ternal d.o.f. were frozen. It was found that the torsion 
of the benzene molecule dominates the chaotic dynamics. 
Without them, there is either anomalous superdiffusion, 
or ballistic motion. 

The total friction at room temperature consists of the 
damping due to the infinite heat bath of the substrate, a 
contribution from the slow dynamics, and one from the 
internal degrees of freedom. These contributions together 
are sufficiently large to account for the high friction found 
in experiments |^ . 

The present results suggest several approaches which 
may make it possible to utilize the internal d.o.f. to affect 
diffusion and friction on a substrate. Molecules which 
are floppy possess many low-frequency vibrational modes 
with strong nonlinear coupling, such as the torsion modes 
in benzene, and consequently, their internal d.o.f. pro- 
duce more noise. Furthermore, as the effective drift and 
diffusion decrease with the time-scale separation param- 
eter, combinations of molecules and substrates where the 
time-scale separation is weaker will produce larger effects. 
Finally, by selectively exciting vibrational frequencies of 
the adsorbate with laser light or by other means, the total 
energy in the internal d.o.f. could be increased, enhanc- 
ing the noise due to the finite heat bath, and therefore 
changing the effective momentum diffusion and friction. 
In this context, large flat photo-active molecules such as 
phthalocyanines are particularly interesting. 
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